Bimodal energy relaxation in quasi-one-dimensional compounds 
with a commensurate modulated ground state 
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We show that the low temperature (T < 0.5 K) time dependent non-exponential energy relaxation 
of quasi-one-dimensional (quasi-lD) compounds strongly differ according to the nature of their mod- 
ulated ground state. For incommensurate ground states, such as in (TMTSF)2PF6 the relaxation 
time distribution is homogeneously shifted to larger time when the duration of the heat input is 
increased, and exhibits in addition a scaling between the width and the position of the peak in the 
relaxation time distribution, it; 2 ~ In (r m ). For a commensurate ground state, as in (TMTTF^PFe, 
the relaxation time spectra show a bimodal character with a weight transfer between well separated 
slow and fast entities. Our interpretation is based on the dynamics of defects in the modulated 
structure, which depend crucially on the degree of commensurability. 
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Non exponential relaxation is the fingerprint of dis- 
order in a system subject to an external perturbation 
0. It is the consequence of a broad distribution of re- 
laxation times. The question of determining the nature 
and spatial extent of correlated objects responsible for 
slow dynamics in spin glasses and glasses, has been de- 
bated for several decades. In spin glasses, the memory 
and chaos experiments give information about exci- 
tations in real space. In glasses there is a consensus 
that structural relaxation occurs on correlated domains 
of a nanometer size 3]. In colloidal suspensions, a di- 
rect three dimensional (3D) imaging of the cooperative 
motion in structural relaxation has been obtained on mi- 
crometer scales 0. A general framework for describing 
spatial and temporal heterogeneities in aging systems has 
been proposed in Ref. 5j. A new class of mesoscopic 
glass, with characteristic scales of micrometers, has been 
recognized among quasi-one-dimensional (quasi ID) sys- 
tems 6j. Heat relaxation at very low temperature T, 
below 0.5 K, in these systems, has been shown to be non 
exponential and to exhibit aging 6] in which the entire 
wide spectrum of relaxation times is shifted to longer 
times with increasing the duration of the very small heat 
perturbation. We suggest here on the basis of heat relax- 
ation experiments that the formation of correlated ob- 
jects in modulated quasi-lD systems, is qualitatively dif- 
ferent in commensurate and incommensurate systems. 

Quasi ID systems exhibit a variety of modulated elec- 
tronic ground states at low-temperature both in the 
charge channel, such as charge density wave (CDW), 
and in the spin channel, such as spin density wave 
(SDW), spin-Peierls (SP) or antiferromagnetic ground 
state. Extensive investigations of the low-T thermody- 
namical properties 0, H @ nave shown that all of 
these systems demonstrate well-known universal features 
of glasses regardless the actual ground state. This in- 



volves an extra contribution to the specific heat C p com- 
ing from low energy excitations, and non exponential 
heat relaxation dynamics for T < 0.5 K (T-range dom- 
inated by these excitations). "Aging" effects are also 
involved, where heat relaxation and correspondingly C p 
depend on the duration of a small heat perturbation (a 
few mK) |E S • Extremely long relaxation times (up 
to 10 4 s) show that these excitations are only weakly 
coupled to the phonon heat-bath [j| , contrary to conven- 
tional glasses. 

The origin of these excitations can be inferred from 
the low frequency dynamics of collective excitations in 
the charge channel. Low frequency dielectric response of 
CDW [lj|, SDW [13 which also has CDW component 
|12j , all having a modulated charge density, is character- 
ized by an overdamped (relaxational) mode slowing down 
with decreasing temperature. This has been recognized 
in the case of CDW and SDW (DW) as the sign of a 
glass transition at the level of electronic superstructure. 
Freezing is naturally generated by screening effects and 
occurs if there are not enough free carriers to screen ef- 
ficiently elastic DW phase deformations. Consequently, 
the characteristic scale of a DW glass is given by the 
phase coherence length, which is of the order of a /xm. 
Below the glass transition temperature T g the low fre- 
quency dynamics and therefore the low energy excita- 
tions are governed by the remaining degrees of freedom - 
topological or plastic deformations such as solitons, do- 
main walls or dislocation loops. 

In this Letter we contrast two isostructural sys- 
tems, (TMTSF)2PF6, which has an incommensurate 
mixed SDW-CDW ground state below 12 K, and 
(TMTTF) 2 PF 6 , which exhibits a charge order transi- 
tion with a ferroelectric behavior at 70 K 0], and a 
spin-Peierls transition at 12 K, both being commensu- 
rate to the underlying lattice (in the following we use 
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FIG. 1: Homogeneous broadening and shift of the heat re- 
laxation with increasing t p (aging) for SDW-PFe. a) Relax- 
ation AT/ATq vs. logi with corresponding stretched expo- 
nential fit exp(-t/r)' 3 , b) S(t) = d(AT/AT )/d(lnt) fitted 
to a Gaussian function G(lnt) with a width w. The direct 
estimate of width at half height shown in b) is 2.17 decades 
in agreement with 1//3 = 2.15 and w = 2.1. 



the notation SDW-PF 6 and SP-PF 6 respectively). Wc 
will show how the degree of commensurability affects the 
very low-T dynamics and discuss these findings in view 
of the properties of solitonic excitations. 

We have measured the low-T heat capacity by stan- 
dard relaxation method. The system is first stabilized 
at a temperature To. Then a small power is fed to the 
sample during a pumping time t p . The temperature 
is increased to the new, slightly higher value T + AT 
(AT/T < 0.05). At t = t p the power is switched off and 
heat relaxation from the sample to the cold source to- 
wards equilibrium at To is recorded as a function of time. 
For all the CDW and SDW systems studied previously, 
below 1 K relaxation starts to deviate from exponential 
on lowering T and also exhibits a strong dependence on 
tp 0- (Fig- Q) ■ A better characterization of non expo- 
nential relaxation can be obtained by the relaxation rate 
S(t) = d(AT/AT )/d(lnt). AT(t)/AT can be written 
in terms of a distribution of relaxation times P(lnr) as 
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where to is the microscopic time. Since T(lnr) varies 
slowly with r we replace the exponential by a step func- 
tion, from what we deduce S(t) = P(bxt) [3. Like in 



FIG. 2: a) An "inverted" Arrhenius plot, 1/T vs. logr m , ob- 
tained for relaxations in thermodynamical equilibrium (r m = 
T eq ) of SDW-PFe. Relaxation rates and their fits are only 
shown for three selected temperatures among five, b) The 
w 2 ~ lnr m scaling is well obeyed as expected from central 
limit theorem. The activation energy is E a = 1.7 K. 



spin glasses there is a "homogeneous" broadening and 
shifting of the relaxation peak with t p (see Fig. HJ-(b)). 
However, since our system eventually reaches thermody- 
namic equilibrium at T + AT for sufficiently long t p pfl , 
aging is then interrupted [l5| . 

Non exponential relaxation can be fitted by a 
"stretched exponential" AT/AT ~ exp [-(t/r)?} 
0, while the relaxation rate S(t) can be fit- 
ted by a Gaussian distribution in hit: G(lni) = 
(1/ViW) exp (— In 2 {t / T m ) / 2w 2 ) , centered on lnr m and 
with a width w. Indeed, the logarithmic derivative of a 
stretched exponential d(exp (— (i/r) /3 ))/d(lnt) is almost 
a Gaussian with r m ~ r and 1/w ~ /3, as shown in 
Fig. ^ The same w can also be extracted directly from 
S(t) curves as the width at half height. 

While the stretched exponential has no known justifi- 
cation other than phenomenological 0, Q , Castaing and 
Souletie argued that a Gaussian distribution in \nt 
can be justified in a framework that associates a hierar- 
chical criterion and a usual scaling approach in the renor- 
malization of distances. Similarly to turbulence they 
use a cascade of events to relate the macroscopic and 
microscopic scales. After iterating the renormalization 
group dynamical scaling takes the form r/r ~ (£/£o) Z j 
where z is the dynamical exponent. The time scale dis- 
tributions are log-normal: P(lnr) = G(lnt) is Gaussian. 



10" 



10° 



10 1 10 2 

t(s) 



10 3 



10 4 



FIG. 3: Bimodal relaxation with redistribution of the spectral 



weight and saturation for t e 



lCT s for the spin-Peierls 



compound SP-PFg. 
figure 



The values of t p are indicated on the 



The variance w 2 scales with the most probable value 
In (t/r m ) 01 following the central limit theorem. It is 
demonstrated in Fig. [21 for SDW-PFg in thermodynam- 
ical equilibrium for three temperatures chosen below 
0.5 K. This yields immediately the temperature depen- 
dence of the stretched exponential exponent f3{T) ~ y/T, 
as (3 ~ 1/w ~ l/VhTr^ (Fig. H>), so that lnr m - 1/T 
indicates a thermally activated process with activation 
energy E a — 1.7 K. At this point we do not intend to 
go deeper into the different models related to (3(T) for 
T —>■ (see ref . @j and references therein) . We stress that 
the w 2 ~ lnr m scaling is quite convincing in the range 
of interest, and also for TaS 3 , a prototype CDW system, 
for which as well as for SDW-PF 6 a glass transition has 
been found at higher T jlOj . 

However, in the case of SP-PFg |9j heat relaxation, 
while still non exponential, exhibits quite different prop- 
erties. Instead of the homogeneous broadening and shift 
of the spectra of SDW-PFg as in a real "aging" ef- 
fects, SP-PF6 demonstrates "discrete bands" of relax- 
ation times. Their distributions are narrower and the 
relaxation times are smaller than in SDW-PF@ (i.e. sat- 
uration of aging is reached within a few tens of minutes 
as compared to days in first system at a similar T). Fig. 
13 shows a bimodal redistribution of the relaxation spec- 
trum indicating that heat perturbation modifies different 
parts of the spectrum in contrast to the homogeneous 
shift of the entire spectrum as in Fig. ^ 

Even in thcrmodynamical equilibrium the S(t) curves 
are not symmetric, indicating at least two different re- 
laxation processes. However the total w (Fig. can 
be obtained, but the variation w 2 ~ lnr m is very pe- 
culiar. We have successfully deconvoluted S(t) in two 
dynamically distinct Gaussian contributions, represented 
by dashed lines in Fig. 0^,. At higher T, the fast contri- 
bution plays the dominant role in the dynamics, but the 




FIG. 4: a) An "inverted" Arrhenius plot, 1/T vs. logr m , 
obtained for the relaxations in thermodynamical equilibrium 
of the spin-Peierls compound SP-PFe, together with a repre- 
sentation of corresponding relaxation rates for slow and fast 
processes and their fits, with E a , s = E a j = 1.0 K b) The 
w 2 ~ lnr,„ scaling does not work for the total relaxation, 
however it might work for each deconvoluted process. 



relative weight of slow modes in the dynamics grows on 
decreasing the temperature. A similar behavior occurs at 
a fixed temperature if we change the duration t p of the 
heat perturbation (Fig. There is an evident weight 
transfer from one to the other subsystem as thermody- 
namic equilibrium is approached. However, without the 
fast entities, there would be almost no response within 
10 2 -10 3 seconds at the lowest T. The distribution of fast 
modes is pushed to 10 3 seconds at 85 mK, and the appar- 
ent activated behavior (Fig. turns into a saturation 
indicating possible tunneling effects (not included in Fig. 
0. 

There are drastic consequences on the proper definition 
and determination of the specific heat C p in these systems 
as it becomes strongly nonlinear ||. From the deconvo- 
luted entities in thermodynamic equilibrium (Fig. 0^), 
it is possible to calculate the corresponding C p , fast and 
slow. Their participation changes progressively from ap- 
proximately 100 % for the fast modes at ~ 650 mK to 
almost 90% for the slow modes below 400 mK. 

It appears that the incommensurate CDW compound 
o-TaS3 shares all its properties related to the shift 
of the dynamics with SDW-PFg. Both are nomi- 
nally incommensurate systems. On the other hand, 
(TMTTF)2Br shares its multimodal dynamics proper- 
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ties with SP-PF 6 Q. (TMTTF) 2 Br being the brother 
compound of SP-PFg is a commensurate antiferromag- 
net. We conclude our Letter by proposing an explanation 
to these disti nctly different dynamics. The models gen- 
eralizing Refs.[h% ll&l20] | are investigated in more details 
inRef. [HI. 

As proposed by Larkin 0] and Ovchinnikov 0], the 
local deformations of the DW in the strong pinning limit 
are bisolitons, obtained by minimizing the elastic and 
pinning energies, and characterized by a ground state 
at energy Eq, separated from another state at energy 
E + AE by a "bounce" state at energy E + AV, there- 
fore defining an effective two-level system. As a qualita- 
tive difference between commensurate and incommensu- 
rate systems plj. one has AE — in the commensurate 
case |22j and AE > in the incommensurate case, so that 
the populations of the effective two level system does not 
couple to temperature variations in the commensurate 
case. We suggest here that, in addition, the degeneracy 
of the classical ground state in the commensurate case is 
lifted by quantum fluctuations that restore a finite heat 
response. The long time dynamics in the incommensu- 
rate case is due to collective effects, as suggested by the 
w 2 ~ In r m scaling on Fig. ^ m agreement with Ref. [2lJ • 
The issue of the long time dynamics in the commensu- 
rate case, involving both quantum and collective effects, 
is mainly an open theoretical question. 

To explain why the spectrum of relaxation time is bi- 
modal in the commensurate case we note that experimen- 
tal data suggest that in the commensurate system there 
exist two entities that relax with approximately the same 
energy barrier but with values of the "microscopic time" 
differing by one order of magnitude. Possible candidates 
for these two entities would be (i) the dipole solitonic 
excitations generated by the strong pinning impurities; 
and (ii) self-induced disorder due to 27r-solitons between 
micro-domains separating the two spin-Peierls ground 
states. Alternatively one may observe that SP-PFg has 
a charge gap Ac — 200 K coexisting with a spin gap 
As ~ 20 K, deduced from the spin-Peierls transition tem- 
perature from the mean field approximation. The coher- 
ence length in the charge channel £c — ^ v f I Ac — 10ao 
is thus much smaller than the coherence length in the 
spin channel £s = Hvf/As — lOOao (where ao is the 
lattice parameter and vf the Fermi velocity). The two 
channels are weakly coupled so that the charge soliton, 
being smaller, relaxes faster than the spin soliton. The 
pinning potential and thus the energy barriers are iden- 
tical in the two sectors so that the activation energies 
are identical but the microscopic times are different, in 
agreement with Fig. 21 

In conclusion we have shown manifestations of corre- 
lated objects interpreted as the collective dynamics of 
a disordered soliton lattice inherent to the modulated 
electronic ground state of DW. We demonstrated that 



a simple scaling feature of the parameters of the relax- 
ation time distribution is well obeyed for the homoge- 
neously broadened response for incommensurate systems. 
By contrast the heat response of commensurate systems 
is bimodal. We stress that our results should be also 
compared to other experiments on the dynamical hetero- 
geneity in the relaxations of supercooled liquids or spin 
glasses |24l Eil Erjf , and to the case of bimodal relaxation 
|23j observed in a fragile glass former at the approach of 
a glass transition. However, the microscopic picture is 
different, in spite of a similar phenomenology. 
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